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We study dark-bright ring solitons in two- component Bose-Einstein condensates. In the Umit of 
large densities of the dark component, we describe the soliton dynamics by means of an equation of 
motion for the ring radius. The presence of the bright, "filling" species is demonstrated to have a 
^3 ■ stabilizing effect on the ring dark soliton. Near the linear limit, we discuss the symmetry-breaking 

Cn ' bifurcations of dark-bright soliton stripes and vortex-bright soliton clusters from the dark-bright 

ring and relate the stabilizing effect of filling to changes in the bifurcation diagram. Finally, we 
show that stabilization by means of a second component is not limited to the radially symmetric 
structures, but can also be observed in a cross-like dark-bright soliton configuration. 

O" 

CN . PACS numbers: 03.75.-Lm, 67.90.+Z 

^ \ Introduction. The dramatic increase of interest in atomic Bose-Einstein condensates (BECs) [l|-[3 ^^^ ^^^ parallel 

(^Jj- significant developments in the field of nonlinear optics [4] have fueled studies of single- and multi-component nonlinear 
Schrodinger (NLS) equations. These operate as a mean- field (so-called Gross-Pitaevskii (GP)) near-zero-temperature 
f^ . model of the condensates in the former case, and as a suitable envelope equation for the electric field propagation in 
Cd ■ the latter. 

In the following we will focus on the realm of multi-component atomic BECs which were realized in a series of 
pioneering experiments p-Q- Among many fundamental effects that have been observed, we highlight longitudinal 
spin waves Q, transitions between triangular and square vortex lattices [9], striated magnetic domains [10|, [ll|, 
phase-separation-induced target patterns [l2[ , tunability of interspecies interactions [l3| and miscibility-immiscibility 
transitions [M]- Furthermore, classes of so-called symbiotic solitary waves (i.e., ones that could not be sustained 
^^ ■ without the nonlinear cross-component interactions) have been theoretically predicted and experimentally observed, 
^ I such as dark-bright (DB) solitary waves [l5l-[l7| and their vortex-bright analog [18]. 

O ■ In this work we demonstrate the existence of a radially symmetric two-dimensional generalization of the DB soliton, 

namely the DB ring soliton whose density profile is shown in Fig. [H In the first, "dark" component the nodal line 
of zero density is accompanied by a phase jump of tt in the order parameter. One-component ring dark solitons have 
been studied e.g. in |19l-[22| and were found to be always unstable towards the formation of vortices. Below, we will 
show that the presence of a second, "bright" component has a stabilizing effect on the symbiotic state. 

Theoretical Analysis. The mean-field model of interest is set up by the two-dimensional (2D) coupled Gross- 
es \ Pitaevskii equations 

CO ■ 

• ; \dt^lJi{x,y,t) = [-|A + V(x,^)+^i|V^i(x,^,t)p+ai2|V^2(2:,^,t)p] V^i(x,^,t) (1) 

O ; idttp2{x,y,t) = [-^A^V{x,y)^g2\Mx^y^t)\^ ^(^i2\Mx^y^t)?]i^2{x,y,t). (2) 

In this dimensionless form, time, length, energy and densities |V^jP, j G {1,2}, are measured in units of cj~^, a^, h^z 
and {2y/27T\ai2\az)~^ , respectively. Here, a;^, a^ denote the oscillator frequency and length in the frozen z-direction, 
respectively (in the x — y plane these are cj^, <^r), whereas an, a22, cli2 refer to the intra- and inter component 
scattering lengths. This scaling leaves us with dimensionless coupling constants gj = ajj/|ai2|, and <Ji2 denotes the 
J^ ■ sign of ai2. Here, we will focus on the experimentally relevant case of ^^Rb hyperfine states {F = l^mp = —1) 
and (F = 2,mF = 1), where all scattering lengths are positive and nearly equal, leading to dimensionless coupHng 
constants of gi = 1.03, g2 = 0.97, cri2 = +1 [l2[- Furthermore, we will assume an isotropic harmonic potential 
V{x^y) = C(j^(x^ + ?/^)/2, and for the numerics fix cj = ujr/ojz = 0.2. As usual, stationary solutions are obtained by 
factorizing il)j{x^y^t) = ex.p{—ifijt)(j)j{x^y), where /ii, /i2 denote the two components' chemical potentials. 

In our first approach, we generalize earlier notions from [15|, |20|, |2J] in order to construct an approximately conserved 
adiabatic invariant of the DB ring, from which a description of its radial dynamics can be inferred. For a highly 
localized DB ring soliton, i.e. for one whose width is much smaller than its radius R^ such an approximate constant 
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of the motion in the regime of large /ii and smah N2 is given by 
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To construct the DB ring energy £^db, we multiply the energy of a one-dimensional DB soliton flSlQ located at distance 
R from the trap center by the circumference 27rR [20|. In this expression, we take into account a scaling correction and 
the rescaled atom numbers Ni^2 = J dxdy\ilji^2{x,y)\'^ . Energy conservation results in the radial equation of motion 
which for a harmonic potential V{R) = uo'^R^ /2 and for low velocities reads 
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The validity of Eq. @ is constrained within radii satisfying jii — V{R) > 0. This equation of motion predicts a stable 
(with respect to radial motions) DB ring equilibrium dit R = Rq (cf. panels (b)-(c) of Fig. [T]) and a breathing radial 
motion around this equilibrium (cf. Fig. [2fa)). 
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FIG. 1. (a) Dark-bright ring soliton's profile at cj = 0.2, /ii = 5, /i2 = 4.07, A^i = 1616, A^2 = 150. (b) DB ring's radius as a 
function of A^i at A^2 = 100 fixed, as calculated from Eq. Q (black) and numerical data (red; gray in the print version), (c) 
DB ring's radius as a function of A^2 at /ii = 5 fixed. 




FIG. 2. (a) Radial oscillations of the DB ring when excited with the breathing mode at /ii = 2, N2 = 100, numerical data 
(red/gray) and harmonic oscillation at the predicted frequency (black), (b) Bogolyubov-de Gennes (BdG) spectrum as a 
function of N2 at /ii = 5 fixed, prediction for the breathing frequency in black. Anomalous modes in the real part shown in 
red/gray, background modes in blue/light gray; purely imaginary modes in red/gray (c) BdG spectrum as a function of Ni at 
A^2 = 20 fixed, exhibiting a stable parameter regime at small Ni. 



^ In Eq. (|3]), we have used the values gi = g2 = cri2 = 1, for mathematical simplicity. We have checked that our numerical results are 
essentially unchanged, if these values are used instead of the experimentally relevant parameters. 



To complement this approach vahd in the regime of large A^i, we also study the DB ring near the linear limit of 
A^i ^ 0, A/'2 ^ 0, where the coupled GPEs reduce to two independent harmonic oscillator Schrodinger equations. 
Stationary solutions of the GPEs have to reduce to harmonic oscillator eigenstates in this limit. To discuss the DB 
ring, we therefore consider linear combinations of the form 

^l)i{x,y,t) =co{t)(po{x,y) ^ ci{t)(pi{x,y) 

for the dark component, where the linear modes cpo = |2,0), cpi = |0,2). Here, \nx^ny) denote the 2D harmonic 
oscillator eigenstates and co,i(t) are complex time-dependent prefactors. This two-mode approximation has been 
shown to be useful in describing one-component ring dark solitons in [23|. For the bright component of the DB ring, 
we numerically find that it reduces to the oscillator ground state |0,0) in the linear limit, and therefore for small 
enough atom numbers we tentatively approximate the order parameter as 1^2 = a/A^|0,0). Note that in this near- 
linear regime the bright component does not yet "fill" the dark ring's zero density region; the DB ring's characteristic 
density profile only emerges in the nonlinear regime of larger atom numbers. Following an approach similar to [25], we 
substitute the above harmonic oscillator ansatze into Eq. ([1]), and use the amplitude-phase decomposition Ck = /3/ce^"^, 
to obtain the consistency condition 
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where a = a^ — ai, and Aqoii = j dxdyipQipi^ ^0122 = /dxd?/(/:?oV^i|'02p etc. denote the nonlinear overlap integrals. 
This algebraic equation for small A^i = /3q + /3f has a single root at /3o = /3i = y^Ni/2, corresponding to the DB 
ring approximately given by ipi oc |0,2) + |2,0), 1^2 oc |0,0) in this regime. Beyond a critical TVf^, two more roots 
emerge, which quickly tend to /3o ^ 0, /3i ^ >/A^ and /3i ^ 0, /3o ^ VNi^ respectively, as Ni is further increased. 
These new solutions correspond to branches of states exhibiting two parallel DB soliton stripes that bifurcate from 
the DB ring branch (cf. Fig. [3]). The effect of this symmetry-breaking bifurcation on the stability of the DB ring and 
its parametric dependence on the bright soliton's atom number N2 will be explored below. 
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FIG. 3. (a) A^i(/ii) bifurcation diagram at A^2 = 20. Solid (dashed) lines denote stable (unstable) branches. (b,c,d) Density 
profiles of the relevant branches (two DB soliton stripes (b), vortex-bright hexagon (c), filled-empty six vortex state (d)) at 
Ni — 100, A^2 = 20. (e) Critical particle numbers at which the two DB soliton branches bifurcate from the DB ring branch as 
a function of A^2: two- mode approach predictions (black) and numerical data (red/gray). The error bars indicate the stepsize 
of the branches' numerical continuation in Ni. 



Numerical Computations. We now turn to a comparison of our theoretical predictions with the numerical results. In 
Figs.[Hb)-(c) the dependence of the DB ring's equilibrium radius Rq is shown as a function of TVi and 7V2, respectively. 
It can be seen that the relevant prediction obtained from Eq. (j4]) is accurate in the proper regime and has the right 
functional dependence on A^i,2- Inaccuracies in the regime of small A^i can be attributed to a breakdown of the 
Thomas- Fermi approximation. For large A^2, we expect repulsive interaction effects between the elementary bright 
solitonic components to become important which are not accounted for in our expression for £^db • The radial breathing 



of the DB ring is shown in Fig. [2]^ a) in fair agreement with the theoretical prediction for the relevant frequency of 
linearization around R = Rq in Eq. (H]). Figures [2fb)-(c) illustrate the dependence on Ni^2 of the frequencies of the 
Bogolyubov-de Gennes (BdG) analysis [2] around a stationary DB ring. When one of these frequencies possesses a 
non- vanishing imaginary part, the corresponding DB rings are dynamically unstable. A key observation here is that 
for fixed A^2, the critical A^i for which imaginary eigenfrequencies arise is larger, the larger N2 is. In fact, for A^2 = 0, 
the dark ring is always unstable [21|, while the presence of the second component yields a stabilization window 
at small A^i; see Fig. EJc). This effect can also be observed in Fig. EJb), showing the BdG frequency spectrum's 
N2 dependence, where a weakening of the dynamical instabilities (i.e., decrease of the corresponding growth rates 
Im(cj)) and stabilization of a number of unstable modes is seen as N2 is increasing. In these BdG plots, the radial 
breathing is represented by the "anomalous mode" of the DB ring whose frequency is in fairly good agreement with 
the corresponding analytical prediction. 

For small atom numbers, we can relate the imaginary modes arising in the DB ring's BdG spectrum, as A^i is 
increased, to symmetry-breaking bifurcations leading to new states. A typical A^i(/ii) bifurcation diagram for a 
relatively small bright atom number of N2 = 20 is shown in Fig. [3l^a). As the dark component's atom number 
increases, branches of two DB soliton stripes bifurcate from the DB ring, destabilizing it. This class of bifurcations 
can be understood within the near-linear framework of the two-mode approximation introduced above. Both in this 
simplified model and in full numerical simulations of the GPE we find that the critical value N^^ where the DB soliton 
stripes' bifurcation occurs increases as a function of A^2, see the very good agreement of these predictions in Fig. 
[3fe). This critical point marks (for small A^2) the termination of the window of stability of the DB ring due to the 
second component referred to above and is associated with the emergence of imaginary modes in its BdG spectrum. 
Subsequent bifurcations from the ring lead to a sequence of polygonal vortex-bright clusters, starting with a hexagonal 
necklace state shown in Fig. ^c). These configurations are natural two-component generalizations of the polygonal 
states discussed in [23|. For small A^2, they are found to be unstable, as they bifurcate from the destabilized DB ring 
branch. The DB soliton stripes also lose their stability to a bifurcating vortex cluster, namely the six vortex state 
shown in Fig. [3fd). Interestingly, in this cluster the bright component only fills four of the six vortices, while two of 
the cores remain empty. For N2 ^ 42 we observe a change in the bifurcation order, with the vortex-bright hexagon 
branch now bifurcating before the DB soliton stripes and thus inheriting the ring's initial stability. 

The above observations show that the presence of the second component parametrically "delays" the instability 
inducing bifurcations in the case of the ring dark soliton. This type of phenomenology is not only present in the 
case of the ring-like state, but can also be seen to persist for other states such as the crossed (diagonal) dark-bright 
solitons of Fig. HJ This can be thought of as a two-component generalization of the always unstable diagonal dark 
solitons of [23|. For this configuration, too, a BdG analysis reveals that for fixed 7V2 > there exists a range of TVi for 
which it is dynamically stable, while beyond the first critical point Ni^, it becomes unstable towards the formation 
of a vortex-bright quadrupole cluster; see Fig. HJc). 
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FIG. 4. (a) BdG spectrum of the crossed DB soliton branch as a function of A^i at A^2 = 30. Same color coding as in Fig. 2, 
oscillatorily unstable modes shown as green/light gray circles, (b) The state's density profile in the stable region (Ni = 10, 
A^2 = 30). (c) An example of a vortex-bright quadrupole cluster for Ni — 70, A^2 = 30. 



Conclusions. In our above exposition, we have illustrated the existence of a series of states in quasi-2D, two- 
component BECs including dark-bright rings, dark-bright multi-stripes, vortex-bright polygons and cross-like dark- 
bright states. For the DB rings and the crossed DB state, we have shown that for a fixed atom number in the bright 
component, the immediate instability of the one-component state is partially inhibited by the presence of the second 
component. The eventual destabilization is caused by the bifurcation of "daughter" states (such as the ones containing 
DB stripes or vortex-bright polygons). For the DB ring it is possible to complement this picture with a (large /ii, 
small A^2) "particle" model which accounts for the existence of an equilibrium ring and its potentially observable 
breathing oscillations. We expect that stabilization by filling should not be limited to 2D realms (see also [17] for 
an example of such stabilization) but could be equally well applicable to one- or three-dimensional contexts. In one 
dimension, a potentially relevant context to consider is the instability of bubbles [26|, while in three-dimensions, the 
study of spherical solitons and of multi- vortex-ring states [3| would be of primary interest. 
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